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Abstract 

The hitherto neglected phonon-exchange interaction between interstitial protons in 
metal lattices is found to be large. It is shown that this effect may give rise to a 
phase of protonic superconductivity, characterized by the formation of Cooper-like 
pairs of protons, in certain metals at high stoichiometric loading. 
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1 Introduction 



The work reported in this Letter was motivated by the observation that, de- 
spite an extensive theoretical literature [1,2,3,4,5,6,7,8] dealing with what are 
often quite subtle effects in hydrogenated metals, no qualitative treatment of 
the phonon-exchange interaction between protons has ever been presented. It 
is hoped that the approximate model outlined here goes some way towards 
remedying this situation. 

The matrix element in the Frohlich model for the (energy-conserving) ex- 
change of a phonon between two charged particles in momentum states {k, —k} 
is (omitting spin indices) 



H t = E \M q \ 2 j= ~ 2 c\ +q c k cl k _ q c- k (1) 



(Ek+q — E k ) 2 — LO, 



.t ~ J 



where u q ~ lob sin(g • a). 

It has been known for a long time [9] that higher-order vertices contribute terms 
of the order of (m/M ion ) 1 / 2 and smaller, so (1), although normally considered 
to apply only to electrons, should also provide a reasonable description of the 
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interaction of interstitial protons in metal lattices. 

At the electron densities prevailing in metals, the screening of a proton is 
well-described by the nonlinear response theory [10,11]. For simplicity, we 
approximate the instantaneous potential due to a proton in a lattice of cell 
dimension a by 

V(r) « (-) ■ e" w where k NL « 2/ a (2) 



Consistency with the classical model of the dielectric effect of the slow moving 
lattice ions convoluted with the adiabatic response of the electron gas requires 

\M I 2 — (3) 



The instantaneous and retarded interactions will be subject to exactly the 
same electronic screening response. The total net interaction between a pair 
of interstitial particles of effective charge Z and mass m is accordingly 



2^2 



Z z e 



, 2 + k 2 



NL 



{q 2 -2k-q)< 



(q 2 — 2k ■ q) 2 — m 2 oj 2 D sin 2 (<f 



C k+q C kC-k-q C ~k 



(4) 



This term is the basis of the BCS theory of low temperature superconductivity. 
There is a qualitative difference from the electronic case however. Whereas for 
electrons of effective mass m* the phonon energies < u q < ujd are typically 
smaller than AE = E k+q — E k everywhere except near the Fermi surface, 
the phonon and quasi-free proton bandwidths are very comparable for a wide 
range of metals because a 2 m p ujD ~ 1- There are thus strong a priori grounds 
for suspecting the existence of paired proton states having binding energies 
considerably greater than those typical for Cooper pairs, in a wide range of 
crystalline environments. The thermodynamic preconditions necessary for the 
actual occurrence of such states will be investigated further in Section 3 below. 



2 Free particle limit 



In order to get a rough idea of the properties that such hypothetical pairs might 
be expected to exhibit, it is necessary to solve (4) explicitly. In anticipation 
of the large binding energy, we neglect the interaction with the lattice. Even 
with this simplification, the two-particle Hamiltonian is still not amenable 
to analytic solution owing to the k-dependence of the interaction potential 
and the periodicity of the phonon dispersion curve. However, if one considers 



2 



only 1=0 states, the two-body system can be transformed into a central force 
problem with spherical Bessel basis set 

^ r) = ^ Cn ^MR ,5) 
n =i rVR 



The resultant k-space matrix can be solved by computer to an accuracy lim- 
ited only by the value of the cut-off momentum N/R. 

With a matrix dimension of N=1000, the calculated energy values and state 
vectors have been found to vary smoothly over a wide range of physically 
realistic values of all the relevant parameters. The dependencies can be sum- 
marized by the following simple generic formulae for the ground state 



E 



-AZe ■ m 



4 • 10" 



kNL \/k NL R 



and (r) 



(6) 



Accordingly, given the assumptions of our model, we conclude that: 

i) the ground state binding energy is of the order of 1 eV and therefore domi- 
nates over the periodic lattice potential, as expected 

ii) the binding energy is proportional to nucleon mass m, (subject torn <C M), 
the bare nucleon charge Ze and the Debye frequency ujd 

iii) there exists a certain critical crystal size of the order of a few nm below 
which no bound state nucleon pairs can form. 

iv) the radial extent of the pair is largely determined by the adiabatic screen- 
ing distance and is therefore orders of magnitude smaller than that typical for 
an electronic Cooper pair. 

Our numerical calculations also revealed excited bound states, provided the 
crystal is large enough. The spectrum is reminiscent of the Rydberg series for 
hydrogen. 



3 Mean-field Hubbard model 



In his analysis of the thermodynamic phases of narrow half-filled bands, Bari 
has shown [12] that (4) admits of a particularly simply approximate form in 
the Wannier representation when a 2 mcj£> is large. Following the approach used 
in that earlier work, we confirm that for a cubic lattice of N interstitial sites 

k 

~ hcjclciCi — —c^jCjCiCi (7) 
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where q is the operator removing a proton at site i, we have omitted the spin 
indices and I\ is the dominant on-site term in (2). In arriving at this grossly 
simplified form, we have used the fact that the major contributions come from 
{k, q} > gT 1 and therefore the {i,j} = terms dominate. 

In view of (6), h ~ leV. The hopping term UjC-Cj in hydrogenated metals 
[7,13] is tiny compared to I\ and h and can be ignored. The effective mean- 
field Hamiltonian for the ith interstitial site is accordingly given by 

H eff = hc\c\ciCi - I 2 {ac\c\ + a*CiCi) - ficjci (8) 



where the off-diagonal order parameter for the superconducting pair state 

a = (c,Cj) must be determined by iteration to self-consistency. 

H is a chemical potential for setting the mean site occupancy (n). 

In view of the identities (c^c^) = (c^c^) = 0, there are just 4 possible Wannier 

site states, with explicit Hamiltonian representation 
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The one-particle partition function, site energies and mean occupancy are 



Z = 2e^ + e- f3E -+e-^ 



2E ± = h-2fi± [(h - 2fx) 2 + 4(/ 2 )>| 2 ] 1/2 



(10) 
(11) 



(n) 



1 

Z 



e^ + 



Eie-P E + E 2 _e-? E - 



E\ + {hfa 2 E 2 _ + (h)*a 2 



(12) 



a : = 



z _ I 2 aE + e-P E + I 2 aE_e-' 3E - 



E\ + (h) 2 a 2 E 2 _ + (I 2 ) 2 a< 



(13) 



Clearly, f3I\ ^> 1 at 3007^. In this "low-T" regime, provided h ~ there is 
a transition from the pure interstitial state characterized by a = to a mixed 
phase with a > above a certain value of (n) < 1. The sharpness of this 
transition increases with \ = hjl\ as shown below: 
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For a more detailed investigation of the phase diagram for this Hubbard model 
and its relationship to the momentum space representation, the reader is re- 
ferred to the work of Kagan and Efremov([14]). 



4 Experimental tests 

The most likely candidates for the proposed phase are the transition metals 
Pd, Ni and Nb. Indeed, some indirect evidence for the onset of protonic su- 
perfluidity in these systems has been reported [15]. In view of the difficulty in 
saturating these hydrides in the bulk, the effect may manifest itself primarily 
at the surface. 

In addition to the normal determinations of superconductivity such as the 
Meissner effect, the exothermy associated with the pairing phase transition 
would be quite considerable and should therefore by readily measurable by 
infra-red or calorimetric techniques. 

It has been pointed out [16] that the proposed proton pairs would possess 
both electric quadrupole (E2) and magnetic dipole (Ml) moments; the corre- 
sponding transitions should be amenable to spectroscopic observation. 
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